Interplay between sublattice and spin symmetry breaking in graphene 
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We study the effect of sublattice symmetry breaking on the electronic, magnetic and transport 
properties of two dimensional graphene as well as zigzag terminated one and zero dimensional 
graphene nanostructures. The systems are described with the Hubbard model within the coUinear 
mean field approximation. We prove that for the non-interacting bipartite lattice with unequal 
number of atoms in each sublattice midgap states still exist in the presence of a staggered on-site 
potential ±A/2 . We compute the phase diagram of both 2D and ID graphene with zigzag edges, at 
half- filling, defined by the normalized interaction strength U /t and A/t, where t is the first neighbor 
hopping. In the case of 2D we find that the system is always insulating and we find the C/c(A) curve 
above which the system goes antiferromagnetic. In ID we find that the system undergoes a phase 
transition from non-magnetic insulator for U < (7c (A) to a phase with ferromagnetic edge order and 
antiferromagnetic inter-edge coupling. The conduction properties of the magnetic phase depend on 
A and can be insulating, conducting and even half- metallic, yet the total magnetic moment in the 
system is zero. We compute the transport properties of a heterojunction with two non-magnetic 
graphene ribbon electrodes connected to a finite length armchair ribbon and we find a strong spin 
filter effect. 
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I. INTRODUCTION 

The most salient electronic properties of graphene and 
its nanostructures are linked to the bipartite nature of 
the honeycomb lattice which is formed by two interpene- 
trating identical triangular sublattices^. It is customary 
to refer to the sublattice as a pseudospin degree of free- 
dom. In this language, the first neighbor hopping is de- 
scribed in terms of a pseudo-spin flip operator, which re- 
sults in the well studied electron-hole symmetric bands in 
graphene, whose wave-functions are sublattice unpolar- 
ized. The pseudospin symmetry becomes a chiral symme- 
try in the continuum limit in which electrons in graphene 
are described with a Dirac Hamiltoniait^ and accounts for 
the lack of backscattering^, the so called chiral tunneling* 
and the absence of an energy gap in two dimensional 
graphene. 

Sublattice symmetry breaking in graphene could 
arise spontaneously, due to some electronic phase 
transition^"—, or due to the coupling of graphene to some 
substrate, like Silicon Carbide^ and Boron NitrideiSrJ^. 
Sublattice symmetry breaking would make it energeti- 
cally favorable for the electrons to stay in one of the 
sublattices, resulting in pseudo-spin order (either spon- 
taneous, or induced). The purpose of this work is to 
understand the interplay between induced pseudo-spin 
order and real spin order in graphene. Magnetic order 
is expected to take place in monohidrogenated graphene 
zigzag edges. Within the standard one-orbital tight- 
binding model of graphene, these edges give rise to a large 
density of states at the Fermi energyi^ which is prone to 
a ferromagnetic inestability when Coulomb repulsion is 



considered within the mean field Hubbard model^''. Den- 
sity functional calculations confirmed the scenario^^^, 
showing that the long-range Coulomb interactions and 
the other atomic orbitals, absent in the Hubbard model, 
do not play a major role in this system. Both the mean 
field Hubbard modeliiii^^^^i^ and DFT calculations 
show that the magnetic phase with zero total spin has a 
gap, which opens due to inter-edge correlations-^-. The 
fabrication of graphene ribbons with ultrasmooth edges 
is now possible by unzipping carbon nanotubes^^"— . In- 
direct evidence of magnetic order in the edges of zigzag 
ribbons is provided by Scanning Tunneling Spectroscopy 
(STS) that can be accounted for within the mean field 
Hubbard modei2£. 

The sublattice degree of freedom plays a central role 
in the magnetic properties of bipartite lattices22r— . Very 
much like an external magnetic field favors one spin orien- 
tation and splits the spin states, an external perturbation 
favors one sublattice with respect to the other and opens 
a gap in the band structure of graphene^ii^. When this 
happens, it is not obvious a priori what happens to the 
edge states, even at the single particle level, and the asso- 
ciated magnetism. DFT calculations indicate that Boron 
Nitride zigzag ribbons with the edge atoms passivated 
with hydrogen are non-magnetio2ii2^ whereas graphene 
ribbons, deposited on Boron Nitride (BN) whose lattice 
parameter is shifted to match graphene , indicate that 
edge magnetism survives^. These DFT calculations sug- 
gest that as the sublattice symmetry breaking potential 
A increases, a phase transition must occur from magnetic 
to non-magnetic edges. Here we address this problem us- 
ing a much simpler description of the electron-electron in- 
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teractions, namely, the mean field approximation for the 
Hubbard model, in the spirit of earlier work for graphene 
with the full sublattice symmetry and on recent work for 
graphene zigzag ribbons without inversion symmetry^. 

The rest of this paper is organized as follows. In sec- 
tion II we present some general theorems regarding the 
properties of the single particle states of the tight-binding 
model for graphene with a staggered potential. In sec- 
tion III we study the interacting model for the case of 
two dimensional graphene and study how the staggered 
potential affects the non-magnetic to antiferromagnetic 
transition. In section IV we study the interplay of mag- 
netic and pseudo-spin order in the case of zigzag ribbons. 
We find that magnetic order and sub lattice symmetry 
breaking give rise to spin polarization of the bands, and 
in some instances we find half-metallic antiferromagnetic 
order. The spin filter properties of this case are stud- 
ied in section V, where wc consider quantum transport 
between two half-metallic zigzag ribbons separated by a 
non-magnetic armchair central region. In section VI we 
summarize our main findings. 



II. SINGLE PARTICLE STATES OF A 
BIPARTITE LATTICE WITH A STAGGERED 
POTENTIAL 



Thus, the chiral symmetry ensures that the spectrum of 
ho has electron hole symmetry. In addition, since ip and 
(j) are eigenvectors of the same Hamiltonian with different 
eigenvalues they must be orthogonal. This leads to: 



or more explicitly 

= J2\^A{^)\ 



(4) 



(5) 



Thus, this lead us directly that eigenstates of ho have 
equal total weight on the two sublattices. In a pseu- 
dospin language, they have a zero expectation value of 
the (Tz pseudospin operator, since the Hamiltonian has 
the pseudomagnetic field (the hopping) in the x, y plane. 

We now turn our attention to the eigenstates of a tight- 
binding hamiltonian with first neighbour hoppings de- 
fined in a bipartite lattice with a sublattice-dependent 
potential which is both homogeneous and traceless, as 
defined by equation ^ . We are going to show that they 
also have electron-hole symmetry. For that matter, we 
represent "Hq in the the subspace defined for a pair of 
eigenstates of /iq, ip and <j) — Uz^P, with energies E and 
—E respectively. We readily obtain 



In this section we consider some quite general proper- 
ties of the single particle states of the Hamiltonian of a 
bipartite sub-lattice with a constant sublattice symmetry 
breaking term: 
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where hAB , h ba and 1 are matrices with dimension given 
by the number of atoms in sublattice A and B. 
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whose eigenvalues are e± = ±y E^ 
sponding eigenvectors v± given by: 



v+ = Cos—tJj + Sin-( 
^ 2 2 

and 
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A. Null sublattice imbalance 

We consider first the case of a bipartite lattice with- 
out sublattice imbalance, so that the number of atoms 
in sublattice A equals those in lattice B: Na = Nb- 
In that case, /iab, hBA and 1 are all matrices of range 
Na = Nb- It can be easily seen that the sub-lattice 
symmetric hamiltonian, ho (or unperturbed hamiltonian) 
anti-conmutcs with the sublattice imbalance operator cr^ : 



V- = Sin—th — Cos—c 
2 2^ 



where CosO = 



(8) 



Thus, if we know (half of) the 



spectrum and the eigenstate of the sublattice symmetric 
problem ho, we can easily build the spectrum and the 
eigenfunctions for the same lattice when a homogeneous 
traceless sublattice Zeeman term is added to the Hamil- 
tonian. 



[ho,(Tz\+ = 



hAB 
hBA 
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B. System with sublattice imbalance 



Since cr^ = 1, it is said that the graphene Hamiltonian 

has a chiral symmetry. As a result, li %jj = \ ) is an 

eigenstate of ho with energy E, we automatically have 
that (j) = aztp is also eigenstate with energy —E 



hip = E^j, h<j>= -Ecj) 



(3) 



The results of the previous section need to be examined 
with care in the special case that E = 0. This certainly 
happens when we consider a system with Na — Nb+Nz, 
where Nz > is a positive integer. In that case the 
dimension of the A and B subspaces is not the same 
and the results of the previous section do not hold in 
general?^. In particular, it has been shown that ho has 



3 




-3 -2-10123 



A/t 

FIG. 1: (Color online), (a) Symbols: energy levels for tri- 
angulene with A'^ = 13 atoms calculated by diagonalization of 
the single-particle model. Lines: energy levels obtained from 
equation (|9]). The degeneracies are indicated in the figure. 



III. ELECTRONIC PROPERTIES OF TWO 
DIMENSIONAL GRAPHENE WITH A 
STAGGERED POTENTIAL 

We now study the interplay between Coulomb repul- 
sion and sublattice symmetry breaking. We model the 
interaction using a Hubbard model in the mean field 
approximation. For symmetric graphene this approxi- 
mation is known to predict a phase transition from the 
non-magnetic gapless state to an antiferromagnetic insu- 
lating state when^ U > Uc = 2.2t. As usual, the mean 
field approximation underestimates the critical U nec- 
essary for the Mott transition. Quantum Monte Carlo 
calculations indicate^ that the transition takes place at 
t/c — 5.3t . In addition, recent work indicates that might 
be a third phase with spin-liquid properties separating 
the non-magnetic state from the magnetically ordered 
phased. In spite of its limitations, the mean field de- 
scription of the Hubbard model can shed some light on 
the possible ordered phases and their electronic proper- 
ties. 



A. Hubbard model and a mean field approximation 



Nz eigenstates ^^:z with E ^ that are are sublattice 

polarized in the majority sublattice, if — ^ ^ ^ ■ This 

are the so called midgap states and play a crucial role in 
the emergence of magnetism in graphene zigzag edges^^ 
and graphene with chemisorbed hydrogen'^^. 

It can be inmediately seen that if ijjz is a zero energy 
eigenstate of h it is also an eigenstate of H=ho + V with 
eigenvalue e = -t--^- Conversely, if ho presents a zero 
energy state sublattice polarized in B, then that state is 
also eigenvector of Ho with energy — y- 

Thus, we can now predict the evolution of the spectrum 
of a given system described by Hq that, at A = has 
midgap states at zero energy as well as pairs of electron- 
hole symmetric states with finite energies ±i?„. The mid 
gap states will split with energy ±A, depending on their 
sublattice polarization, and the finite energy states will 
evolve as 



e±{A) = ±^El + —. (9) 

In order to illustrate this result, we have computed the 
single particle spectrum of a triangulene^^ with Na — 
7 and Nb = 6 atoms. The evolution of the spectrum 
as a function of A is shown in figure ([T]). We compare 
the result of the numerical diagonalization with those 
extrapolated from the spectrum of ho and, expectedly, 
find perfect agreement. 

This calculation shows that, in structures with a larger 
number of midgap states, the midgap shell will remain 
half- full (when counting the spin), and interactions are 
expected to favor large spin configurations. 



The extended Hubbard model reads: 

i,s 

nifUii (10) 

i i 

where cj^ creates an electron in atomic site i with spin 
s =t,i, i' stand for the first neighbors of i, Tz{i) — +1 
if i belongs to the A sublattice and —1 otherwise. We 
only consider the half-filling case, where the number of 
electrons equals the number of sites in the lattice. For 
a given filling, the ground state properties of the model 
depend on two dimensionless parameters A/t and U/t. 
We explore the properties for a spin-coUinear mean field 
approximation, where the U term is approximated by: 

Vmf = +U ^ rii^iriii) + {ni^)nii (11) 

i 

where (riis) is the average of the occupation operator of 
site i with spin s, calculated with the many-body ground 
state of the mean field Hamiltonian: 

(n,,) =^/„(a|n,,|a) (12) 

a 

where — 0,1 is the occupation of the single particle 
states I a) that diagonalize calculated with the ground 
state of the mean field Hamiltonian Ho + Vmf- Since 
the potential Vmf depends on the eigenstates of H + 
Vmf, both the potential and the eigenstates need to be 
computed self-consistently. We do this by iteration. 



4 



B. Mean field approximation for 2D graphene 

We now describe the electronic properties of Hubbard 
model for the two dimensional honeycomb lattice with 
a staggered potential within the mean field approxima- 
tion. In this case, we can take a minimal unit cell with 
2 atoms, A and B and assume that the mean field in all 
unit cells is identical, which permits to use Bloch theo- 
rem to represent the mean field Hamiltonian in the basis 
set A't,Bt,Ai,B i: 



H 



i?i 



where each element is a 2 by 2 matrix 



2 i U{nAi) 



and 



and 



f + U{nA^) 

rik) 



fik) 



fik) = t(l + 
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(13) 



(14) 



(15) 



(16) 



accounts for the first neighbour particle hopping. In our 
numerical determination of the self consistent occupa- 
tions {riis) we have taken an unit cell of 4 atoms. We 
have verified that our mean field solutions in this ex- 
tended unit cell do not present inter-cell modulations of 
the charge density. 

We have explored the phase diagram defined by U/t 
and A/t and we find 3 types of solution, shown in figure 
@: 

1. For U < Uc{A) the system is non magnetic and, 
except for A = 0, a band insulator. The case of 
A = and U < Uc is the well studied paramagnetic 
semimetal phase. 

2. For U > Uc{A) the system is an antifcrromag- 
netic insulator. The lack of inversion symmetry 
produced by the sublattice symmetry breaking re- 
sults in a splitting of the spin bands, in contrast 
with the standard A — case. This is shown in 
figure © 

3. For U — Uc{A) the system is a half-semimetallic 
antiferromagnet. For one spin channel the system 
is insulating and for the other is semimetallic. 

It is apparent that, as A increases, the critical Uc in- 
creases. Expectedly, the magnetic order has to overcome 
the single-particle gap opened by the staggered poten- 
tial. Interestingly, the mean field approximation describe 
a magnetic transition between two insulating states, the 
non-magnetic insulator and the antiferromagnetic insu- 
lating phase, which can be interpreted as an excitonic 
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FIG. 2: (Color online). Phase diagrams for 2D-graphene 
(left) and a zigzag graphene nanoribbon with = 48 atoms 
in the unit cell (right) with stagger potential (A) using a 
mean field Hubbard model at half filling. The dark region 
with Amit{U) > A{U) > Ac{U) correspond to the spin half- 
metallic phase in the graphene ribbon. In the case of 2D- 
graphene, this region is reduced into a single critical line sep- 
arating non-magnetic and antiferromagnetic insulating states. 
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FIG. 3: (Color online). Band structure of 2D-graphene with 
stagger potential (A) in the mean field extended Hubbard 
approximation for U > ?7c(A) 



insulator transition. Given the large values of Uc{A)/t 
this ordered electronic phase is not expected in graphene. 
The predictions of this theory should be tested in cold 
atomic gases confined in optical lattices^ or in artificially 
paterned honeycomb lattices in two dimensional electron 
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IV. ELECTRONIC PROPERTIES OF 
GRAPHENE ZIGZAG RIBBONS WITH A 
STAGGERED POTENTIAL 

We now study the case of zigzag graphene ribbons for 
which we find that magnetic order could happen at low 
values of U/t even for finite A. The width of the ribbons 
is characterized by N, the number of atoms in the unit 
cell. Importantly, one of the edges is formed with A 
atoms only, the other being made of B atoms only. Thus, 



pseudospin polarization implies charge accumulation in 
one edge and depletion in the other, ie, the formation of 
an electric dipole. 

The electronic structure of graphene ribbons has been 
widely studied in the A = limit, both for the U = Oi^iii 
and the finite U casesi^"— i^Sii^. The most prominent fea- 
ture of their electronic structure is given by the flat bands 
associated to edge states. At A = J7 = 0, these bands are 
located at the Fermi energy, giving rise to a large density 
of states at the Fermi energy. Not surprisingly, Coulomb 
repulsion results in a magnetic inestabilityi^ correspond- 
ing to the formation of magnetic moments in both edges 
while the bulk-atoms remain almost spin unpolarized. It 
turns out that the inter-edge spin correlations are antifer- 
romagnetic, as expected from the Lieb theorem^''. Thus, 
for a given spin orientation, there is charge accumulation 
in one of the edges and charge depletion in the opposite. 
This results in a spin-resolved pseudo-spin polarization, 
or spin-dipole^^. Here we are interested in the interplay 
between pseudo-spin polarization, driven by the A term 
in the Hamiltonian, and the spin polarization, which en- 
tails a spin-resolved pseudo-spin polarization, driven by 
the Coulomb repulsion U. 



A. Non interacting bands 
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FIG. 4: (Color online). Electronic structure of zigzag ribbons 
with [/ = , for A = and A = 0.2t (left and right panels) 
for 2 different ribbon widths, A = 40 (top) and, TV = 80 (bot- 
tom) . For the sake of clarity, we only plot the 2 higher energy 
valence bands and the 2 lowest energy conduction bands. 



We first review the effect of the staggered potential on 
the non-interacting bands, studied by Qiao et al^. At 
A = two almost flat bands, associated to edge states, lie 
at the Fermi energy. As A becomes finite (and positive) , 
the bands at B edge are red-shifted and those at B edge 
are blue-shifted, resulting in a band-gap opening. This is 
seen in figureQ for two ribbons with A^ = 40 and A^ = 80 
atoms, for A = (left columns) and A = 0.2t (right 
column). We also notice the low energy bands are quite 
similar for the A^ = 40 and A^ = 80 ribbons, whereas the 
gap between the higher energy bands is reduced for the 
wider ribbon. This is consistent with the fact that lowest 
energy bands are edge states, relatively insensitive to the 
width of the ribbon, in contrast with higher energy bands 
made of quantum confined bulk states^. 

Thus, for finite A and U — Q graphene zigzag rib- 
bons are band insulators with pseudospin polarization 
that features two fiat bands corresponding to the highest 
occupied and lowest un-occupied bands. For A > 0, the 
bands corresponding to both t and J, spins in the B edge 
are occupied, whereas those in edge A are empty. As we 
show now, these bands are prone to magnetic inestability, 
not-unlike in the case with A = 0. 



B. Effect of Coulomb repulsion 

We study now the interplay between pseudospin po- 
larization and Coulomb repulsion. For that matter, we 
use again the mean field approximation for the Hubbard 
model, as described in previous wor k-'^'^'-'^'^d^-^^i^^-'^'^ . Nu- 



merically found solutions present magnetization at both 
edges equal in magnitude and opposite in sign. Thus, 
there are two equivalent ground states: niA — —ms > 
and ruA = ^ruB < 0. The corresponding energy bands 
for the ribbon with A^ = 48, and A = and U = t 
are shown in figure ([IK). The magnetic order results in 
a band-gap opening. The spin f and \. bands are de- 
generate. The magnetic moment at the edge atoms is 
m = ±0.13. The charge per atom is the same all over 
the unit cell, 1 electron per atom. 

When the sublattice-symmetry breaking potential is 
finite and below a critical value Ac(C/), we still find mag- 
netic order in the edges with antiferromagnetic coupling, 
and zero total moment, even if the charge is no longer 
the same for both edges. The electronic properties of the 
magnetic ribbon with pseudo-spin polarization are differ- 
ent on several counts. First, the bands are spin split, as 
a natural consequence of the lack of both time-reversal 
and inversion symmetries. The evolution of the energy 
bands, as we increase A, is shown in figure © for the 
A^ = 48 ribbon with U = t. This figure can be under- 
stood as follows. The solution has — —ms > so 
that the t band is occupied (empty) in the A (B) edge. 
Conversely, the J, band is occupied (empty) in the B (A) 
edge. As A is turned on, the B bands are red-shifted and 
A bands are blue shifted. For the t bands, this implies 
that the band gap closes, since the valence A bands move 
upwards and the conduction B bands move downwards. 
Conversely, the gap opens in the J, channel. 

As shown in figure ^ and discussed below, as A in- 




FIG. 5: (Color online). Lowest energy bands for for zigzag 
ribbon with A'^ = 48 atoms in unit cell and U = t, for different 
values of A: (a) A = 0, (b) A = O.OSt, (c) A = 0.2t (d) 
A = 0.3t. Only the 2 highest energy occupied bands and the 
2 lowest energy empty bands, per spin channel, are shown. 
Blue (red) stands for t (V) bands. 



creases the magnetic moment at the edges are depleted 
and, eventually disappear when A > Ac(U). Remark- 
ably, the gap in the t channel closes for Amit(C^) < 
Ac{U), yet the gap is finite in the I channel. Thus, 
the combination of pseudo-spin polarization and antifer- 
romagnetic order makes the system a half-metallic anti- 
ferromagnet in a region of the A, U phase space (see dark- 
violet region in figure ([fright panel). Notice that this is 
different from the ferromagnetic half-metallic phase pre- 
dicted for graphene ribbons in the presence of a trans- 
verse electric fieldii, for which the total magnetic mo- 
ment is different from zero. 

Finally, we note that the gap of the non-magnetic case 
in figure ([5|)d for A > Ac is significantly smaller than 
A. This is due to the renormalization of the bands due 
to Coulomb repulsion. Basically, the occupied bands are 
blue shifted with respect to the empty bands, reducing 
the size of the gap. 



C. Phase Diagram 

In figure ^ we show the phase diagram defined by 
U/t and A/t for a ribbon with = 48 atoms, calculated 
within the mean field approximation at half-filling. Ear- 
lier work^i has addressed the phase diagram defined by A 
and the electron density. The diagram in figure ^ has 2 
phases regarding the magnetic order: non-magnetic, for 
A > Ac{U) and antiferromagnetic otherwise. The re- 
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FIG. 6: (Color online). Properties of = 48 ribbon with 
U — t as sl function of A. Top panel: edge magnetization. 
Middle panel: edge charges. Bottom panel: gap for spin f 
and 4- 



suits are very similar for ribbons with different widths. 
In contrast with the 2D case, for A = the critical U 
for the edge is zero. This makes zigzag graphene rib- 
bons suitable systems for the observation of magnetism 
in graphene and the possible effect of sublattice sym- 
metry breaking more relevant. Expectedly, the critical 
Ac{U) is an increasing function of U, or in other words: 
the larger the single particle gap, the strongest the in- 
teraction U required to drive the magnetic inestability. 
Spin polarization requires promotion of electrons across 
the single particle gap, from the occupied B to the empty 
A edge, ie, the formation of a magnetic exciton conden- 
sate. The difference with the A = case stands on the 
size of the the single particle gap, which is vanishingly 
small (but not zero) in finite width ribbonai^. Interest- 
ingly, the magnetic exciton condensation scenario already 
takes place in the case apparently conducting A = U = 
ribbon, when U is turned on. 

For a fixed value of A, as the value of U is increased the 
system undergoes a phase transition from a non-magnetic 
band insulator, to a half-metallic antiferromagnet and 
then to an insulating antiferromagnet. The fact that 
interactions can drive the system from insulating from 
metallic is quite exotic and differs from the usual Mott in- 
sulator scenario, in which interactions drive a band metal 
insulating. In the phase diagram we mark the insulator 
to metal transition when the smallest energy gap, at a 
given spin channel, is 100 times smaller than the gap at 
A = 0. As shown in figure ^ br U = t, the gap at 
A = is 0.07t. Thus, for that particular value of /7, we 
declare the system conducting when the gap is 7 10~'^t. 
Variations upon this criteria yield quantitative changes in 
the metal to insulator transition line, but it is always the 
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case that Amit{U) runs along and below the magnetic 
phase transition line A(.{U). 



V. GRAPHENE ZIGZAG RIBBONS WITH 
STAGGERED POTENTIAL AS AN IDEAL SPIN 
INJECTOR 

Interestingly, the predicted conducting phase for 
Amit{U) < a < Ac{U) is a half-metallic antiferromag- 
net. In this section we study the spin transport proper- 
ties of a tunnel junction where the electrodes are made 
of such half-metallic antiferromagnets and the barrier is 
made of semiconducting armchair graphene ribbon. 

In figure© we consider three possible situations, all 
of them with A = 0.25i , U = t. The first and second 
cases feature two antiferromagnetic electrodes with mu- 
tually parallel and antiparallel magnetizations, respec- 
tively. Based on the band structures of these infinite rib- 
bons, shown in figure, we expect the tunnel conductance 
to be completely depleted when the magnetic moments of 
the different electrodes are anti-parallel. In the last case 
we consider two ferromagnetic electrodes. The bands of 
the infinite ribbon with ferromagnetic coupling between 
the magnetic edges reflect the conducting and spin unpo- 
larized character of the system at the Fermi energy (set 
at E=0 eV in the three cases studied). 

The conductance of the system is calculated within 
the Landauer formalism and the Green's funtion method 
in a system consisting of an armchair semiconducting 
graphene nanoribbon with an on-site Coulomb potential 
connected to two spin-polarized zigzag nanoribbons with 
both stagger and on-site Coulomb potentials, as those 
studied in the previous section. Except for the atoms 
at the interface with the zigzag ribbons, the Hubbard 
U term is not able to spin-polarized the amrchair cen- 
tral region, as shown in the figure ([Tj. The short length 
of the tunneling barrier justifies also the neglect of spin 
relaxation which are expected in longer samples. 

To compute the transmission function T{E) along the 
armchair nanoribbon we adopt a partitioning method 
as implemented in the ALACANT(Ant.U)^ transport 
package, where the system is divided into three parts22i2i, 
namely, the central (C) part which consists on the arm- 
chair ribbon connected to a small part of the electrodes, 
and the left and right semi-infinite staggered zigzag 
nanoribbons (L and R). 

The transmission probability can then be obtained 
from the Caroli expression!^ 

T^{E) = TT[gliE)TR{E)gciE)TL{E)]^ (17) 

where GcaiE) = [zl - Hc^a - ^rA^) - ^lA^T^ is 
the Green's function of the central region which contains 
all the information concerning the electronic structure of 
the semi-infinite leads through the self-energies {TiR^cr{E) 

and T^lAE)), and TrAE) = A^lAE) - ^^^(i;)], 

^L,cr{E) = i[Y.R^c!{E) — e]^^(£')] are the coupling ma- 



trices containing the information about the coupling of 
the central region to the leads. 

In Fig. [51[a-f), we have computed the spin-polarized 
conductance for the three systems shown in Fig. [71 The 
top and middle panels correspond to the cases with and 
without stagger potential respectively. The bottom panel 
of Fig. [5] shows the conductance polarization P = G^^ — 
Gj. /G-i- -I- Gj. X 100 at different energies with (violet line) 
and without (green line) stagger potential. 

The results obtained for the spin conductance can be 
easily inferred by looking at the band structures of the 
three regions in Fig. [71 The two antiferromagnetic cases 
without stagger show a large gap in the conductance 
due to the semiconducting behavior of the three regions. 
For the ferromagnetic case, there is a finite conductance 
near the Fermi energy due to the metallic nature of the 
electrodes. In this case, the evanescent modes coming 
from both electrodes penetrate into the central region 
and overlap due to the short length of the armchair rib- 
bon. 

When both electrodes are antiferromagnetic with mu- 
tually parallel magnetization, the system transforms into 
a spin-half metal where both valence and conduction 
bands show the same spin polarization. In this case, the 
polarization of the conductance is the same at both sides 
of the Fermi energy. The case featuring antiferromag- 
netic leads with mutually antiparallel magnetization also 
shows a spin half-metallic phase in both leads but with 
opposite spin polarization around the Fermi level. Thus 
the spin polarized current injected from one electrode is 
always reflected by the opposite electrode resulting in a 
zero conductance around the Fermi energy. If both elec- 
trodes are coupled ferromagnetically the system becomes 
metallic and the valence and conduction bands cross at 
the Fermi level but with a different spin polarization. 
This makes the conductance polarization to change sign 
at each side of the Fermi level. 



VI. SUMMARY AND CONCLUSIONS 

Our main findings can be summarized as follows: 

1. We have shown that for a bipartite lattice with a 
site-independent pseudo-spin Zeeman term the re- 
sulting spectrum has electron hole symmetry, ex- 
cept for mid-gap states that arise in lattices with 
a different number of sites in the two sublattices. 
We have shown that, if the A = Hamiltonian has 
an eigenstate (j) with energy E, it has also an eigen- 
state '0 with energy —E and the Hamiltonian with 
finite A has two eigenstates, linear combination of 

$ and if with energies ± \J ^ + E^ 

2. If (f) localized in the A (B) sublattice is an eigenstate 
with energy E ^ oi the A = Hamiltonian, 
then it is also an eigenstate of the Hamiltonian with 
finite A and energy A (—A) 
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FIG. 7: (Color online) Band structure of the leads and cen- 
tral region for the three cases studied: (a) parallel antiferro- 
magnetic electrodes, (b) antiparallel antiferromagnetic elec- 
trodes and (c) parallel ferromagnetic electrodes,. The cen- 
tral region is the same for the three cases studied, namely, a 
semiconducting armchair graphene nanoribbon with a local 
Coulomb potential {U ~ t). 
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FIG. 8: (Color online) Spin conductance through a finite 
armchair ribbon connected to two spin-polarized zigzag rib- 
bons with (a-c) and without (d-f) stagger potential. The bot- 
tom panel shows the spin conductance polarization for each 
magnetic ordering at the electrodes in the presence (violet) 
and absence (green) of stagger potential. 



3. Two dimensional graphene with finite A undergoes 
a transition to an antiferromagnetic state with spin- 
split bands, for U > J7c(A). The critical C/c is an 
increasing function of A. At the transition between 
the non-magnetic insulating state {U < C/c(A)) and 
the magnetic state U > Uc, the system is a half 
metallic antiferromagnet. 

4. Zigzag graphene ribbons with a finite sublat- 
tice symmetry-breaking potential below a criti- 
cal value (A(f7) < Ac{U)) can undergo a transi- 
tion from non-magnetic insulators to spin-polarized 
half-metallic antiferromagnets in the presence of 
Coulomb interaction. The fact that interactions 
can drive the system from insulating to metallic 
differs from the usual Mott insulator scenario in 
which interactions drive a band metal insulating. 

5. Zigzag graphene ribbons with stagger potential A 
slightly below the Ac{U) are predicted to be ideal 
spin injectors. The spin transport calculations car- 
ried out in this work show high spin polarization of 
the conductance (w 100%) around the Fermi level 
when used as spin injectors in a tunnel junction. 
This indicates that, if the suitable substrate that 
yields the right A is found, graphene zigzag rib- 
bons could act as half-metallic spin injectors. 

This work has been financially supported by MEC- 
Spain (FIS- and CONSOLIDER CSD2007-0010). We ac- 
knowledge fruitful conversations with Jeil Jung. 
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